RKKY interaction in a disordered two-dimensional electron gas 
with Rashba and Dresselhaus spin-orbit couplings 
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We study theoretically the statistical properties of the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction between localized magnetic moments in a disordered two-dimensional electron gas with 
both Rashba and Dresselhaus spin-orbit couplings. Averaging over disorder, the static spin sus- 
ceptibility tensor is evaluated diagrammatically in the mesoscopic (phase-coherent) regime. The 
disorder-averaged susceptibility leads to a twisted exchange interaction suppressed exponentially 
with distance, whereas the second-order correlations, which determine the fluctuations (variance) 
of the RKKY energy, decay with the same power-law as in the clean case. We obtain analytic 
expressions in the limits of large/small spin orbit interactions and for equal Rashba and Dresselhaus 
couplings. Beside these limiting cases, we study numerically the variance of the RKKY interaction 
in the presence of pure Rashba spin-orbit coupling. Our results are relevant for magnetic impurities 
or nuclear moments embedded in III-V two-dimensional heterostructures or in contact with surface 
states of metals and metal alloys, which can display a sizable Rashba spin-orbit coupling. 

PACS numbers: 75.75.-c, 71.10.Ca, 75.70.Tj, 71.23.-k 



I. INTRODUCTION 

The Ruderman-Kittel-Kasuya-Yosida (RKKY) 
magnetic interaction between nuclei 1 or magnetic 
impurities^ in bulk metals is determined by the non- 
local spin susceptibility of the conduction electrons. 
A given localized magnetic moment in contact with 
a clean electron gas at low temperature produces an 
electronic spin polarization which is oscillatory (with 
period half of the Fermi wavelength) and decays in 
magnitude as a power law of the distance from the 
impurity. Other distant magnetic moments interact 
with such polarization pattern. In the presence of 
weak non-magnetic disorder for the electron gas, the 
configuration average of such long range interaction is 
exponentially suppressed beyond the mean-free path. 4 
This is due to the randomization of the oscillatory tail of 
the interaction but does not imply that the interaction 
strength is exponentially suppressed. The magnetic 
interaction is still long-range and the variance is found 
to decay with the same power law of the clean case.^* 6 - 

Similar results are also valid for electrons confined in 
two dimensions^— Furthermore, in semiconductor het- 
erostructures, the long range nature of the RKKY in- 
teraction is enhanced in the presence of nonanalytic 
corrections to the spin susceptibility* 1 ^— and a mag- 
netic ordering transition of the nuclear spin system was 
predicted ] 10 ' 11 In the following, we neglect these interac- 
tion effects (becoming relevant for low electron densities) 
and focus instead on the role of the spin-orbit coupling, 
which is responsible of anisotropic magnetic interactions. 
Spin-orbit scattering by nonmagnetic impurities is im- 
portant to determine the anisotropy field of bulk spin 
glasses^ and can be treated by standard diagrammatic 
techniques.— iZriLiii Its effect on the statistical properties 
of the magnetic interactions in two dimensions was stud- 
ied in Refs.[8|andl9|, which in particular obtain a crossover 



from Heisenberg to Ising behavior. 

In this paper we assume that spin-orbit scattering from 
the impurities is not present, but we consider trans- 
lationally invariant spin-orbit interactions due to the 
band structured The effect of the Rahsba spin-orbit 
coupling^ was examined in Refs. [l8l - |2ll and a simple 
form of a twisted RKKY exchange interaction was de- 
rived at large distance.— A discussion including also 
the Dresselhaus spin-orbit interaction^ can be found in 
Ref. |23|. All such studies neglect the effect of disorder but 
this assumption is not always justified, since the spin- 
orbit length can be larger than the mean-free path in 
realistic conditions. It is the main purpose of this pa- 
per to examine in detail the interplay of disorder with 
Rashba and Dresselhaus spin-orbit couplings. 

Besides being applicable to III-V heterostructures, 
where the electron gas mediates interactions between nu- 
clear spins, our results are relevant for isolated mag- 
netic impurities on metal surfaces. In recent experi- 
ments, the RKKY interaction of a pair of such adatoms 
could be directly probed as a function of their relative 



position 
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On the other hand, surface states with 



a sizable Rashba spin-orbit coupling exist, as first ob- 
served on Au(lll)i 26 ' 27 Much larger spin-orbit split- 
tings are found in metal alloy s 28 ! 29 and several meth- 
ods exist to modify the strength of the spin-orbit cou- 
pling in these systems i^2r— Finally, Rashba and Dres- 
selhaus spin-orbit coupling could influence the prop- 
erties of low-dimensional magnetic semiconductors as 
(In,Mn)As and (Ga,Mn)As where the RKKY interaction 
is mediated by holes £2r£l Indeed, a large spin-orbit in- 
teraction is present in two-dimensional heterostructures 
for holeS ) 16 ' 38 i 39 but the dominant spin-orbit coupling 
is in this case cubic in momentum instead of being 
linear . 16 ' 38 ' 40 Although our results are not directly ap- 
plicable, it would be possible to treat this relevant case 
in a similar fashion. 



The paper is organized as follows: In Sec. |TT]we discuss 
the main definitions and the general formalism of our 
work. The diagrammatic expressions for the disorder- 
averaged quantities of interest are introduced. An ex- 
plicit calculation of the susceptibility tensor can be found 
in Sec. IIIII and the second order-correlations of the sus- 
ceptibility tensor are obtained in Sec IIVI for some impor- 
tant limiting cases. More generally, the variance of the 
RKKY interaction has to be evaluated numerically and 
the calculation is presented with many details in Sec. [V] 
for pure Rashba spin-orbit coupling. Finally, Sec. lVII con- 
tains our conclusions and a number of technical points are 
discussed in Appendices |A"1Id1 



II. GENERAL FORMALISM 

We consider a system of localized moments (e.g. mag- 
netic impurities or nuclear spins) interacting with a free 
electron gas. Restricting ourselves to a specific pair of 
such impurities located at positions Ri, R2 with spin 
operators Ii, I2 we have 



H = Y J H el (i) + jY, S(Ri)-Ii, 



(1) 



where H e i(i) is the single-particle electron Hamiltonian 
and the second term in Eq. ([J} is a contact interaction 
with the two magnetic impurities, expressed in terms of 
the electron spin density (in units of Ti/2) 



S(R) = J2 6 ( r * - R ) <*i 



(2) 



with Ti, <Xj the position and Pauli spin operators of elec- 
tron i. We consider here point-like magnetic impurities, 
i.e., with spatial extent much smaller than the Fermi 
wavelength^ We also assume an isotropic interaction of 
the Heisenberg type, e.g., describing the hyperfine inter- 
action of conduction electrons in III-V semiconductors, 
but the following discussion can also be easily adapted to 
an Ising-like interaction, more appropriate for holes42 

In a coordinate system with x, y respectively along the 
crystal axes [100] and [010], the single-particle electron 
Hamiltonian has the following form 

P 2 

H e i = Va(py(j x -p x (jy) + f3(p x o x ~pyay) + V(r), (3) 

which includes the Rashba and Dresselhaus spin-orbit 
couplings and the disorder potential V^(r). Equation ([3]) 
has a simpler form if x, y are along [110] and [110], as 
assumed through the rest of the paper 



H e i = h a-Pya x - a + p x ay + V(r), (4) 

where a± = a ± /3 define the following spin-orbit lengths 



1 



Notice that Eqs. (|3]) and ([J} have the same form if /3 — 
0, since the Rashba spin-orbit interaction is rotationally 
invariant in the xy-plane. In this case, a single spin-orbit 
length A = l/m|a| appears. 

By making use of second-order perturbation theory in 
J, the effective RKKY interaction between Ii and I 2 is 
obtained as follow*^ 

H 12 = - J 2 Y.^ihjXij (^1,^2), (6) 

in terms of the static spin susceptibility tensor of the 
electron system Xij(Ri,R2), given by^ 

X«(Ri,Ra) =-j-J ([^(Ri,t),^(R 2 )]-)e- r ' t dt, 

(7) 
where (. . .) denotes a thermal average, [a, b]± = ab ± ba, 

77 = 0+, and S(R() = e _ * i: i ir ' ,(i) *S(R)e* E i H « ,W3 *. 
As in Refs. [a and [3 it is convenient to use the Matsub- 
ara technique and by analytic continuation to imaginary 
frequencies Eq. ([7]) is rewritten 

(Ri,Ra) - -r^Tr[a l G w „(R 1 ,R 2 )a J G w „(R 2 ,R 1 )]. 



Xij 



(8) 
In Eq. ©, T is the temperature, uj n = (2n + 1)ttT are 
Matsubara frequencies (we set ks — % — 1 in the fol- 
lowing), the trace is taken over the spin indexes, and 
G w (Ri,R2) is the single particle electron Green's func- 
tion 

^(Ri.Ra) = (Ri|(iw - H el + 6 F )- 1 |R 2 ), (9) 

where ep is the Fermi energy and |R) are eigenstates of 
the position operator. 

Since the disorder potential V(r) is unknown, we per- 
form an average over its possible realizations (denoted 
by an overbar) . We introduce the following notations for 
the susceptibility tensor 



Xw(R) =Xu(Ri,R2), 
and its second-order correlations 



XijXi'j'(R) = Xy(Ri,R2)x l 'j'( R i 1 R 2 



(10) 



(11) 



where we could set Ri = and R 2 — R = R cos ip e x + 
Rsunpe y , since translational invariance is restored. 
Here, et are unit vectors along the coordinate axes. 

B y assuming the simplest case of (5-correlated disor- 
der V{r)V(r') = {mry^ir - r'), Eqs. (PS5J> and (TTT]) 
can be calculated to leading order in 1/kpl, where 
kp = mvp = \J2raep is the Fermi momentum and 
I = vpT is the mean-free path. The standard diagram- 
matic technique can be applied by including the spin- 
orbit coup ling in the d isorder-averaged Green's function 
G W (R) = G w (Ri,R 2 ). In the self-consistent Born ap- 
proximation G u (R) is given by 



m a± 



(5) 



G W (R) 






G„(k)e 



ik-R 



(12) 





44-46 



FIG. 1. Diffuson (left) and cooperon (right) diagrams con- 
tributing to the second-order correlations of the RKKY inter- 
action tensor given in Eq. (|15[) . 



A± 3> t, Ao(q, u) is given byss 
X D (q.,u) = 

-2m(. 2 

\ u u vnor_ \ 

— ia+q x —ia_q y m(a^_ + a 2 _) 
which allows to evaluate Eqs. (fTo]) and (ITS1) . 




(18) 



when 



,44 46 



III. ASYMPTOTIC EXPRESSION OF THE 
SUSCEPTIBILITY TENSOR 



GJk) 



i (u> + 22^) - is_ + £f + a_k y (T x - a + k x a y 



[i(u 



It 



k 2 

2m 



£f] -a-^-a+fc? 



(13) 

We refer to Sec. IIIII for an explicit evaluation of Eq. (TT2")) 
in the asymptotic limit R 3> 1/kp. 



A. General formulas for Xsj and XijXi'j' 

The leading contribution to Eq. (fit))) has a simple form: 
Xy(R) ~ -r^TV^G^fR^CJ-R)], (14) 

n 

which corresponds to a single empty bubble diagram. 
This expression applies in the diffusive limit, when R 3> I 
and vertex corrections can be neglected. On the other 
hand, XijXi'j'C^) is obtained to leading order as a 
sum of two series of diagrams, involving either diffuson 
or cooperon propagators (see Fig. [T^iLS The equality 
of diffuson and cooperon diagrams for vanishing exter- 
nal magnetic field is known for other fluctuation cal- 
culations (in particular, for the universal conductance 
fluctuations^) and we show in Appendix [Al that this oc- 
curs also here. Therefore, we can write the final result in 
terms of the diffuson propagator Dr(w) as follows 

/ \4 />oo 

x^aW(R) = ^/ wdu J2 ^rH^r'H 

IT JO , i 

The diffuson Dr(oj) is a 4 x 4 matrix given by 
1 r dq 



D R (u) 



2mr J (2tt) 2 



-***■[!- X D {q,w)]-\ (16) 



with 



A^(q, W ) 



rlk 



where the indices (A, v assume values 0, x, y, z (eg is the 
identity and <J x ,y,z are the Pauli matrices). In the diffu- 
sive limit q£, \uj\t <C 1 and for small spin-orbit couplings 



We first study the disorder-averaged susceptibility ten- 
sor, given by Eq. (|14[) . Clearly, xij(R) is immediately ob- 
tained from the real space representation of the Green's 
function [see Eqs. (|T2"1) and (fT3|) ] and we start from the 
angular integral 



G W (R) 



kdk 
(2^ 



rf^G^kje-^™*'-*', 



(19) 

which can be performed with the method of the station- 
ary phase if kR 3> 1. For small spin-orbit couplings, 
this condition is realized when R 3> 1/kp because the 
largest contribution to the integral is from the region of 
the Fermi surfaces. A sum of two terms is obtained, cor- 
responding to ip p ~ if and ip v ~ f + Tt ((pis the azimuthal 
angle of R) 



G " (R) "?i Wf- 



2nk 
~R 



,^i(kR-ir/4) 



[i(u 



2t J 



~ £; + £f} 2 - k 2 al 



, (20) 



where a v = a_ sin fe x — a+ cos if e y . The main con- 
tribution to the fc-integrals is from the poles located 
(to leading order in the spin-orbit coupling) at k± — 
kp ± ma v + i (u + 2S1H) j VF , If w > 0, only the term of 
the sum proportional to e+'C^-^/ 4 ) gives a nonvanishing 
contribution 



G r [i ,>o(R) 



y/2-KkpR 



E 



e »(fc±K-7r/4) 1 

2 



1±^- 



o. 



(21) 
The physical meaning of this result is rather transpar- 
ent since Re[fc±] are the Fermi momenta of the two spin 
branches along the direction of R. Furthermore, Eq. (f2"Tj) 
involves the spin projectors on the eigenstates of a v ■ cr, 
determined by the spin-orbit coupling of Eq. Q when k 
is along R. By combining this result with the expression 
at uj < we find 



(17) G w (R)~-i 



, (sgn a;) TO i( sgnw ) (k F R-n/4) -{\u\ + l/2r)R/v F 



y/2nk F R 
cos[ma v R) + i- 



■ sin(TOa ¥ ,i?) 



(22) 



which is in agreement with the asymptotic form of the 
exact Green's function with only Rashba spin-orbit cou- 
pling (see Ref. LL8|), except that qp — y/Kp + m 2 a 2 is re- 
placed here by kp. This difference arises because in the 
derivation we neglected corrections of order a\, which 
implies that our result is only accurate for R <C kp A±. 
However, if the condition kp\± ^> 1 is valid (which is 
often satisfied in practice), Eq. (f2"2"j) becomes inaccurate 
only on a length scale much larger than min{A+, A_}. 

By using Eqs. (fl4|) and (|22|) . the following result for 
the susceptibility tensor is obtained at T — » 

^ ^y ^,^^ (23) 

wheieTZfj D — hTi[aiU(TjU'] and U appears in the second 
line of Eq. J25J 



C7 



ARmOLta-tT 



(24) 



Clearly, [/ is a rotation operator of angle 2ma v R around 
dtp. This gives that 7Z RD is a rotation matrix, which 
generalizes the twisted exchange interaction of Ref. [T^. 

The effect of the impurity potential in Eq. (|2"5|) is the 
factor e~ R / 1 . Since the disorder-averaged calculation is 
valid in the diffusive limit R 3> I, the final result for xJJ 
is exponentially small and it is necessary to consider the 
second-order correlations XijXi'j' ( see following sections). 
However, Eq. (|23|) is still useful in the ballistic case, when 
the expression for Xij (R-) is formally identical to Eq. (fT4|) 
if the r —¥ oo limit is taken in the Green's functions. 
Therefore, Eq. (|2"5]l is also valid for R <C £. In this case, 
both disorder-averaging and the exponential suppression 
can be dropped. 

Finally, we note that the averaged charge response 
Xoo(R-) can also be simply calculated from Eq. (f2"2"j) and 
is the same as in the absence of spin-orbit coupling. 
This is at variance with recent results of Ref. 47J, where 
anisotropic Friedel oscillations were obtained, which can 
also display beatings as function of R. We attribute this 
discrepancy to our approximation of small spin-orbit cou- 
pling. The effects of Ref. |47J are quadratic in a± at small 
spin-orbit couplings and thus are only relevant for dis- 
tances much larger than A±, when Eq. (|22|l becomes in- 
accurate. 



As discussed in Appendix [Cl the only dependence of 
rjij^'j' is from the ratios R/X± and the direction of R. 
This is similar to the results of Refs. [5J and 0, where 
the ratio R/L SX} , appears instead of R/X± (L s . . is the 
spin-orbit diffusion length, in the presence of spin-orbit 
scattering) . While it is possible to evaluate explicitly the 
tensor rjij^'j' as function of R/X± (see Sec. |V] for pure 
Rashba spin-orbit coupling), we start discussing Eq. (f2"5)l 
in some limiting cases that can be solved analytically and 
have a direct physical interpretation. 



A. Small spin-orbit coupling: £ <C R <C A± 

In this case, since the relevant momenta and frequen- 
cies of Xp>(q, u>) are of order q ~ 1/R and lot ~ (£/R) 2 , 
we can neglect the second line of Eq. (fT8|) (proportional 
to the spin-orbit couplings) and the diffuson (fll))) can be 
evaluated analytically as£ 



Dg(*) 



*o(f 



IixwitI' 1 



(26) 



where Kq(z) is a modified Bessel function of the second 
kind4£ Finally, Eq. (JTSJ) gives 



XijXi'j>(R) = 3n4R4 8ij8i'j> 



(27) 



The variance of the RKKY energy immediately follows 
from XijXi'j' [see Eq. ©] 



{Hl2)2 -3^ {Il ' h)2 



(28) 



Therefore, the interaction between the magnetic mo- 
ments is of the Heisenberg type. This is appropriate for a 
system without spin-orbit coupling^ where the nonlocal 
polarization created by the impurity at the origin is ev- 
erywhere parallel to Ii. Equation (|2"5)l reveals that such 
polarization is long-ranged and decays in magnitude as 
1/R 2 . However, its strength is randomized in sign and 
the average susceptibility is exponentially suppressed, see 
Eq. (E2D. 



IV. SPECIAL LIMITS OF THE 
SECOND-ORDER CORRELATIONS 



The general structure of XijXi'j' 
Eq. (fl"5)) can be represented as 



following from 



XijXi'j' {"■) 



2m 2 
3vr 4 i? 4 



f]ij,i'j' 



(25) 



where t]ij%'ji is defined in Eq. (JC2J). We find that such 
a tensor is of order unity, which implies that the RKKY 
interaction decays in magnitude as 1/R 2 . 



B. Large spin-orbit coupling: I <C A± -C R 

In the opposite limit of a large distance between the 
magnetic impurities (relative to the spin-orbit precession 
lengths A±), the spin-orbit contribution to Eq. (TT8l) is 
much larger than the first line. Therefore, [1 — Xp>]oo is 
much smaller than all the other matrix elements. Since 
the diffuson is obtained from (1 — Xjj)^ 1 we can approx- 
imate Eq. (fl"8t as 



DgM 



2irmTp 



■ SfiOOvO, 



(29) 



which gives 



XijXi>j'(~&) 



m x c 



(30) 



The variance of the RKKY energy reads in this case 



H 



12 J 



J m 
6tt 4 .R 4 



t2t2 

1 1 1 2 



(31) 



Therefore, the magnetic interaction decays with the same 
power law 1/R 2 as in the absence of spin-orbit couplings, 
see Eq. (|28|) . However, the variance becomes independent 
of the orientations of the two magnetic moments since I 2 
are constants. This fact becomes clear taking into ac- 
count the spin relaxation of the electron spins induced 
by the spin-orbit interaction, which occurs on a length 
scale X± . The direction of the electron polarization in- 
duced by the first impurity is parallel to Ii close to its 
location (i.e., at the origin) but is fully randomized at a 
distance R S> A±. Since the spin polarization induced 
in R by the first impurity is finite but distributed on 
average with spherical symmetry, Eq. pip can be sim- 
ply understood. Notice also that Eq. (|3Tj) is reduced by 
a factor 1/4 with respect to the maximum variance of 
Eq. (EH]). 



C. The case \a\ = \/3\ 

We consider here Rashba and Dresselhaus spin-orbit 
couplings with equal strength, when a particular com- 
ponent of the spin is conserved.— For definiteness, we 
assume a — j3 (i.e., a_ = 0) such that a y is conserved, 
while a similar argument is valid if a + = 0. This situation 
is most simply treated using an exact unitary transforma- 
tion since the electron Hamiltonian © can be expressed 
as 



Hei = c/ f hr- + y ( r ) ) U + 

2vn ' 



(32) 



independently of the disorder potential V(r), where U — 
e -i{ma + )xc y ig tne p era tor of Eq. ([23]). Therefore, for a 
given disorder potential, xy(0, R) is obtained from the 
result without spin-orbit coupling after a rotation of I2 
(and a small shift of the chemical potential cf)- Such 
rotation is around e y by an angle 2ma+Rcos<p. 

Since this exact property holds for every disorder real- 
ization, it still holds after averaging. The result for the 
susceptibility tensor is in agreement with Eq. (|23p while 
the second-order correlations are 



XijXi'j'CR) 



2m 2 
3tt 4 .R 4 



qjRDqjRD 



(33) 



with IZ-if the same rotation matrix of Eq. (|23[) . evalu- 
ated at a_ = 0. The variance of the energy is 



(#i 



2J 4 m 2 

37T 4 i? 4 



(Ii-ft 



RD-, 



(34) 



The same result can be obtained by direct evaluation of 
Eq. (fi"5|) . see Appendix [5] 

If the spin orbit couplings are not exactly equal, 
Eq. (|33j) is only valid for R <C A_ and, by increasing R, 
a smooth crossover to Eq. (f30|) is realized when R ;g> A_. 
Notice however that our treatment is only valid if R is 
within the phase-coherence length of the electron gas. 

The regime close to a = /3 was also examined in 
Refs. and Ho, where it was found that several lin- 
ear response functions have a non-analytic behavior. For 
example, the electrically induced spin-polarization of a 
phase coherent sample of size L is zero or finite depend- 
ing on the order of the limits L — } 00 and a — > j3^ 
Similarly, we find in our case that XijXVf nas two differ- 



ent forms, Eqs. (|30p and (f3"3"|). depending on the order of 

the R — > 00 and A_ — > 00 limits. 



V. VARIANCE OF THE RKKY INTERACTION 
WITH ONLY RASHBA COUPLING 

In the general case, we could not obtain analytic results 
for Eq. (1151) . We then resort to direct numerical integra- 
tion and consider here the specific case of pure Rashba 
spin-orbit coupling (j3 = 0). Since the Rashba spin-orbit 
coupling has rotational symmetry in the xy-plane, we as- 
sume in the following R = Re x . Furthermore, results for 
pure Dresselhaus spin-orbit coupling can be simply ob- 
tained applying a spin rotation. The general form of the 
correlation tensor XijXi'j' is given in Eq. (|25|) in terms 
of the tensor Tjij^if , which uniquely depends on the ratio 
R/X (since for Rashba spin-orbit coupling A± = A). De- 
tails on the evaluation of the functions f]ij^ y are given 
in Appendix [Cl while we find it more transparent to dis- 
cuss here the variance of the RKKY interaction energy 
instead of the full correlation tensor. We first discuss 
specific contributions to the variance, corresponding to 
definite orientations of I1.2 in the classical limit. A gen- 
erally applicable expression Eq. (|5T|) is obtained at the 
end of this section. 



A. Ii and I2 parallel to e y 

As a first example, we consider the contribution to 
the variance involving only I\ y and I^y-, i.e., due to the 
XyyXyy component of the correlations tensor: 



(Hl2)'< 



y-y 



2J 4 m 2 
3tt 4 .R 4 



^lyliy 



(35) 



From Eq. (J25J), it is clear that we defined Aq = f]yy.y y , 
which is plotted in Fig. [2] and is a monotonically decreas- 
ing function of R/X, from 1 to 1/4. Therefore, it correctly 
interpolates between the two limits of Sees. IIV Al and 
IIVBI We note that the asymptotic value 1/4 is reached 
with good approximation already at moderate values of 
R/X > 2. The reason is that the asymptotic value is 




.-._. 1(T<5A 



-o 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

R/A 

FIG. 2. (Color online.) Main plot: dependence of Aq of 
Eq. (|35|) on R/X (red solid curve). The horizontal dashed line 
is the asymptotic value 1/4. The coefficient A' of Eq. (|45p 
almost perfectly coincides with Ao- The blue dashed curve is 
the small difference SAo — Ao — A' (multiplied by a factor of 
100). The inset shows Ao — 1/4 (dots) on a semilogarithmic 
plot. A linear fit gives ln(A - 1/4) ~ 1.73 - 2.57(R/X). 








R/A 



FIG. 3. (Color online.) Main plot: the four solid red over- 
lapping curves show the dependence on R/X of the rotation 
angles $, $', and $ ± defined in Eqs. ([38]>. (|44)l . and (|49)l . 
respectively. The dashed line indicates the linear dependence 
<!> = R/X — maR. Inset: the small differences $' — $ (purple, 
solid), $ + — $ (blue, long dashed), and $ _ — <I> (red, short 
dashed) are plotted as functions of R/X. 



approached exponentially, as it can be seen in the inset 
of Fig. [21 A similar behavior was found for spin-orbit 
scattering, where the exponentially small corrections in 
the limit R/L so , 3> 1 were explicitly evaluatedj££ 



B. Ii along e y and 1 2 in the perpendicular plane 

As a second case, we consider terms involving only I\ v 
and the x, z components of I2: 

AJ in 2 / T 2 t2 

= 3-7T 4 /? 4 ly^y x 'V x '■ 2x "1 Vyz,vzJ-2z 

+Vyx,yz[hx,hz] + ), (36) 



{Hl2Y 



where [a,b]± = ab ± ba. We do not provide here ex- 
plicit expressions for the components of r]ij,i>j> enter- 
ing Eq. (|36[) but they can be found as discussed in 
Appendix [C] In particular, they follow directly from 
Eq. (|C2p . Equation (fMf can be directly applied to clas- 
sical magnetic moments if Ii along e y and I2 in the xz- 
plane, and we can simply set [hx, l2z]+ = 2l2xl2z hi this 
case. More generally, this is not justified since Ii 2 are 
angular momentum operators and Eq. (|36[) represents a 
specific contribution to the total variance Eq. (|51"j) . 

Eq. (|5rj)) can be simplified with a rotation of I2 along 
e^ by defining 



S2X 

S 2 z 



cos $ — sin $ 
sin $ cos $ 



/ 2 , 



where $ satisfies the condition 

^Tjyx.yz 



tan 2$ 



Vyz,yz T]yx,yz 



(37) 



(38) 



The plot of $ as function of R/X is given in Fig.[3J After 
this transformation, the anticommutator term drops out 




FIG. 4. (Color online.) Dependence on R/X of A x (purple, 
solid), A y (blue, long dashed), and A z (red, short dashed). 



and Eq. (|5fj)l becomes 



(#12) S 



2J 4 m 2 

37T 4 i? 4 



\AzI\yJ2x 



AxllSl), (39) 



where the coefficients A^, and A z smoothly interpolate 
from to 1/4 and are plotted in Fig. @] 

If the magnetic moments Ii and I2 are interchanged, 
we can proceed in a similar way. From Appendix [D] we 
have that r) X y,zy — —Vyx.yzi implying that the final ex- 
pression is simplified by rotating Ii by an angle — $. This 
gives 



2J 4 m 2 

o 4d4 (A x S lz I 2y + A z S lx I 2 y)- (40) 



{H12Y 



Since Ii and I2 are transformed with opposite angles, the 
two rotated reference frames differ by 2$. From Fig.[3]we 
see that 2$ ~ 2maR, which is reminiscent of the twisted 
exchange interaction of Ref . [l8| and Eq. (I2li)) . 



C. Ii and I2 perpendicular to e y 

Consider now the contributions from XijXi'f if a H the 
indexes are different from y. As it turns out, not all such 
components are independent. Besides the symmetries 
discussed in Appendix [D] we find rj X x,xx = t]zz,zz and 
Vxx.xz = Vzz,xz- The three other relevant components of 
the tensor are T} XZtXZ , r) xz , zx , and T] X x,z Z - Therefore, the 
following expression is obtained 



(tfu) 



J 4 m 2 

37 



+(r?x a ,x :c +r?x 2 ,x 2 )(Ii 2 x + A 2 J(4 2 x + 4) 

+ (»?xx, Z ^-»?x Z)Z x)[-flx.-fl z ]-[i2x,^]-}, (41) 



where we defined 

J 2 -I 2 

IX %Z 

[-l-ix-! -*izj-f- 



K: = 



A* 



^]xx,xx Vxz,xz 
^Vxx,x 



^Vxx,xz 
'}xx,zz 1 Vxz 



(42) 



(43) 



and V^ is the transpose of Vi. By a suitable rotation 
of I 2 as in Eq. (p?Tf and of Ii by an opposite angle, the 
expression V^ MV2 can be diagonalized. However, the 
angle is not exactly $ as before, but a slightly different 
value <&' which satisfies 



tan 4$' = 



4rfe 



1]xz,xz T]xx,xx TJxx.zz T]xz,zx 



(44) 



As seen in Fig. |3l the plot of $ and $' are practically 
indistinguishable. The dependence of the small difference 
5$' = <&' — $ on R/X is shown in the inset of Fig. 03 In 
terms of the rotated magnetic moments S^ ~ S j, Eq. (T4"Tj) 
takes the form 

_ 2J 4 m 2 
x 2 ~~ 3vr 4 i? 4 ' 



(#i 2 ) s 



7 , By Pl x j <->lzJ±P2xi ^2zJ± 
± 
, 4/ /q/2 c/2 , c/2 c/2\ 

+MStS'l + S? z Sg)). (45) 

The coefficients A y , 2?+, and B~ are plotted in Figs. SJ 
[3 and IH respectively. Finally, the value of A' is very 
similar to A : its plot and the small difference SAq = 
A — A' are shown Fig. [5J As discussed in more detail in 
Appendix [C] the fact that the angles $ and $' and the 
amplitudes Aq and A' are very similar can be ascribed to 
the same reason, i.e., the diffuson matrix Z)r(w) can be 
diagonalized at each value of R/X with a rotation which 
is almost independent of the frequency u). The rotation 
angle is approximately equal to $. 



D. General case 

We consider now an expression valid for generic orien- 
tations of the magnetic moments. Such generic formula 
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FIG. 5. (Color online.) Dependence on R/X of B~£ (purple, 
solid), By (blue, long dashed), and Bt (red, short dashed). 



BT 




FIG. 6. (Color online.) Dependence on R/X of B x (purple, 
solid), By (blue, long dashed), and B~ (red, short dashed). 



for the variance contains the sum of Eqs. Q35p. (|59")l . (|40l) . 
and (|45)l . but an additional term A(Hi2) 2 appears, in- 
volving the remaining six independent elements of the 
correlation tensor rjij^/ji. This last contribution is ex- 
pressed in compact form as 



a(h 12 ) 2 



where we defined 



J A m 2 tt^^t 



37T 4 i? 4 ^ 



T7- _ [ V-iyi hx\± 
*l± - I [r t ] 

\*-iy ) *iz\ ± 



V 1 1 ± M±V 2 ±, (46) 



(47) 



U [ Vxx.yy ± Vxy,yx Vxz.yy ± Vxy.yz \ /^g\ 

\ Vxz,yy T Vxy.yz 1]zz,yy ^ T^zy^yz J 

As before, the matrices A^i can be diagonalized with 
opposite rotations of Ii and I2. The angles are different 
for the two terms entering Eq. (|46p . They are given by 



tan2$ d 



^\J]xz,yy ^ T]xy,yz) 



Vxx,yy ^ T]xy,yx ~r f]zz,yy ^ Vzy,yz 



(49) 



and are plotted in Fig. [3J As seen, $ ± are both almost 
identical to $ of Eq. ([35]). The differences (5$ ± = $ ± - $ 
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R/A 



FIG. 7. (Color online.) Thick black solid line: plot of 
Bo = 2B^ as function of R/X. This can be compared to 
the coefficients 2Bf (purple solid lines, see also Fig. [SJ and 
2B~ (short dashed red lines, see also Fig. [6]). Thick black 
dashed line: plot of Ao — Bo, which can be compared to the 
Ai coefficients (long dashed blue lines, see also Fig. [4|. 



the expression, one can set B+ ~ B+ = Bq/2 and 
A x>y>z ~ An — Bq. The error introduced by these approx- 
imations is seen in Fig. [Jj where we plot Bo and An — Bo 
together with the corresponding coefficients 2Bf and Ai 
of the previous Figs. [5] and [U Equation (|5T1) then assumes 
the following simple form 

gj^ [5o(Sx • S 2 ) 2 + (A - B )SJSi] , 

(53) 
at R/X < 



(H12) 2 



which explicitly interpolates between Eq. 
1 (when A = Bo = 1) and Eq. fl3TJ at R/X > 1 (when 
Aq = 1/4 and Bo = 0). Equation (|53")) can also be used to 
approximate the quantum case, if both the commutator 
and anticommutator coefficients are set equal to Bo/2. 
This can be justified from the fact that all the 2B~f co- 
efficients decrease from 1 to on a similar scale of R/X. 
However, the approximation Br ~ Bo/2 is significantly 
worse than for Bf (see Fig. [7J . 



are plotted in the inset of Fig. [3l which shows that 5$ 
are in general very small and only increase in magni- 
tude at large values of R/X. However, as discussed in 
Sec. HVBi the angular dependence of the magnetic mo- 
ments becomes unimportant in this limit. In terms of 
the rotated magnetic moments S i ~ Si, we can rewrite 
Eq. |gB|) as follows 



A(ff 12 ) 2 = 



2J 4 m 



\B X l^lyi ^lz\±[^2yT ^2z\± 

+Bf[Sf x , St v ]±[S£,S±] ± }, (50) 



37T 4 i? 4 ^ 



where the coefficients B+ z , and B~ z are plotted in Figs. 
[5]and[6l respectively. 

Finally, we can summarize the final result for the vari- 
ance as follows 

2J 4 m 2 r 

(-H12) — o 4j->4 "I AqS 1x S 2x + A x (S ly S 2z + S lz S 2y ) 

+ J2 B t[Siy,Su}±[S 2y ,S 2z } ± + c.p.}, (51) 



where c.p. indicates the remaining terms, obtained by 
cyclic permutations of the indexes xyz. Two such terms 
are of the form AoS^S^i (i — y,z) and the other six 
correspond to the coefficients A y , A z and By,Bf. Notice 
that in Eq. (|5~1T) we neglected the small difference between 
A and A' as well as between Si and S^, S { . If we do not 
use this approximation and we apply Eq. (fBTj) to spin- 1/2 
magnetic impurities we obtain 



(H12) 2 



37T 4 i? 4 



An 



A' a + J2 ( A i- 2 KSuS^ 



t—x,y,z 



(52) 
by using S' ly = S iy = I iy . 

For the opposite case of classical variables, the com- 
mutator terms in Eq. (1511) vanish. To further simplify 



VI. CONCLUSION 

We discussed in this paper the influence of non- 
magnetic disorder and spin-orbit couplings on the RKKY 
interaction in two dimensions. Our analysis focused on 
the effect of the Rashba and Dresselhaus spin-orbit in- 
teractions, rather than spin-orbit scattering^ We exam- 
ined both the ballistic and disorder-averaged suscepti- 
bility. Furthermore, we characterized the second-order 
correlations of the susceptibility tensor. 

For the susceptibility tensor, we extended the result 
found in Ref. [l8| and obtained that the RKKY interac- 
tion has the form of a twisted Heisenberg exchange for 
generic couplings a, /3. An interesting question is if such 
form of twisted exchange coupling is still valid at large 
spin-orbit couplings or distances, beyond the order of ap- 
proximation considered here. In this regime, anisotropic 
Friedel oscillations and the presence of beatings were ob- 
tained in the density response^ 

At distance larger than the mean-free-path, the 
disorder-averaged RKKY interaction is exponentially 
suppressed. Therefore, it has to be characterized through 
the second-order correlations of the susceptibility ten- 
sor, which decay with the same power law 1/i? 4 of the 
clean case. We obtained the tensor structure explicitly in 
a number of interesting limits and presented a detailed 
study with pure Rashba spin-orbit coupling. In this case, 
the final result Eq. (|51l) is rather involved but the main 
features are well reproduced by a simpler approximate 
formula, Eq. (|5"3")) . We obtain that: (i) also for R >• £ 
the RKKY interaction (in this case, its variance) can be 
expressed in terms of rotated magnetic moments. The 
relative angle deviates from the linear dependence of the 
clean case (see Fig. [3]). (ii) The regime of large spin-orbit 
coupling is realized at moderate distance. As exempli- 
fied in Fig. [21 the coefficients of Eq. (|51l) approach their 
asymptotic values already at R/X > 2 and the corrections 
are exponentially small. 



We provide next a summary of relevant length scales 
and experimental conditions to which our results ap- 
ply. First of all, it is required for all the distances in- 
volved in the problem to be smaller than the phase- 
coherence length L v (which can be of the order of a few 
micrometers^), the size L of the sample, and the thermal 
diffusion length 6 - (L T = ^D/ttT, with D = v F r/2 the 
diffusion constant). We also generally assumed through- 
out the paper that 1/fcp is the smallest length. For mag- 
netic semiconductors this implies that, since R 3> l/k F , 
the free carriers should be provided mainly form external 
doping and not by the magnetic ions. The impurity cal- 
culations are generally performed in the diffusive regime, 
which implies R ^> £. Finally, effects of the spin-orbit 
coupling are relevant when the spin-orbit lengths become 
comparable to the distance R. Therefore, we can sum- 
marize all these conditions as follows: 



{L v , L, L T } > R, A± > t > 1/fcp. 



(54) 



Additionally, our result for xij(R-) can be applied to the 
ballistic case. This is because the Green's functions ap- 
pearing in the bubble diagram Eq. (fT4]) become equal to 
the unperturbed ones in the limit r — > oo. Therefore, 
Eq. ([25)1 becomes a valid expression for \ij (R) (without 
disorder-averaging) if 



mm{i p , L, L T , £, \±k F } > R, A± > l/k F . 



(55) 



In conclusion, our study is relevant to characterize 
magnetic structures and spin-glass properties in two- 
dimensional systems. However, instead of having sev- 
eral magnetic impurities, a single isolated magnetic mo- 
ment could be considered. The spin polarization pat- 
tern induced in the electron gas could be revealed for 
surface states by scanning tunneling microscopy (STM) 
techniques! 24 ' 25 STM studies also exist in magnetic 
semiconductors^ but optical methods might represent a 
more viable option. 36 In particular, spatially resolved op- 
tical imaging of the spin polarization is possible in semi- 
conductor quantum wells^ Optical pumping could also 
be exploited to polarize the central magnetic impurity 
or a small region of nuclear spins in a chosen direction. 
Finally, an alternative approach might be provided by 
spin-grating techniques^ 3 - to obtain the spin susceptibil- 
ity in momentum space. 



ACKNOWLEDGMENTS 



We acknowledge useful discussions with O. Chalaev, 
M. Duckheim, D. L. Maslov, and M. J. Schmidt. This 
work was supported by the Swiss NSF and the NCCR 
Nanoscience Basel. 



-9^- 



I i i 



■<- 



a' 
a 



->- 



b' 

b 



I I I 



->■ 



b' 



FIG. 8. Diffuson (top) and cooperon (bottom) ladder dia- 
grams corresponding to Eqs. (|A1|) and (|A2[) . respectively. 



Appendix A: Cooperon and diffuson contributions 

to XijXi'f 

We discuss here in more detail the derivation of 
Eq. (TTol) for XijXi'j'- The final result involves both diffu- 
son and cooperon diagrams (see Fig. [IJ , constructed from 
ladders schematically illustrated in Fig. [3] The diffuson 
ladder corresponds to the analytic expression 

— ~ / II tS [G* (ki) • • ■ G Ul (k„)] a6 



(2tt) 2 (mr)"" 1 ./ 11 (2tt) 2 

x [G U2 (k„ - q) . . . G U2 (ki - q)] 6 , , ,(Al) 
while for the cooperon ladder we have 

(2^) 2 (mr)»-ii 2(^2 [G Ul (ki)...G Ul Qs n )U 

x [G W2 (q - ki) . . . G U2 (q - k n )] a , b , .(A2) 

In Eqs. (jAip and (IA2[) the Green's functions in the 
square parentheses are 2x2 matrices (in the spin space), 
a, b, a', b' are spin indexes, and (n — 1) is the number of 
impurity lines (n > 2). We now make repeated use use 
of the identity 4 ^ 



1 

Sabred — -Z 2 J°>]qc[°>]c 

for the diffuson ladder and 

3 

2 



1 3 

Sabred — - 2_J[° '^]ac[o fj,)o 



(A3) 



(A4) 



n=o 



for the cooperon, where a^ — a y <J^. This allows to sim- 
plify ([ATI) and (1A"2J) as follows 



— Yy^ aa '\. av \ b ' b j j2^i' 



^L--*»- R [(jr D (q,w)) n U(A5) 



rriT 



^2Wl]aa'W v ] 



b'b 



f/q 



iq-R, 



■ w [(A c (q, W ))"]^,(A6) 



where Xo(q, ui) is defined in Eq. ([TT]) and 
X ^ (q ' " } = 2^7 / (2^ Tr[ ^ G - (k) " iG - (q " k) 
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In the above equations, ui = \lji — W2I and G^(k) is the 
transposed Green's function. 

By summing Eqs. (|A5j) and (|A6[) over n we ob- 
tain (mT) 2 [a fl ]aa'[o'u]b'bD^(uj) for the diffuson and 
(mr) 2 [aj l ] a a'[& I j]b'bC^l' (uj) for the cooperon case. The 
cooperon propagator Gr(w) is defined similarly to 
Eq. ([To]) , with Ac instead of Xo- The full series of fluc- 
tuations diagram is than easily obtained 



XijXi'j' (&) = Yl 12 ( 



mT) 4 T 2 



x [L>^(w)D^ (w)Tr(o'icr A1 cri'cr y ')Tr(o' : ,-(T A1 'cr J '0' y ) 

+C^(a;)C^'( W )Tr( ( 7^tcTj^)Tr(cT^>^,)]. (A7) 

We can finally take the limit of zero temperature. This 
amounts to replace T 2 Y2 Ul u 2 ~^ 2^ /o°° ^duj, using that 
the integrand depends on uj = |wi — W2I and is zero if io\ 
and W2 have the same sign. Furthermore, we can use 
(TycrJ <jy = —ui (for i = x,y, z) in the two cooperon spin 
traces and obtain 



XijXi'j'(R)= 



(mr) 



2tt 2 



LJ 



duj 










C&»C^(«)] 



x 



Tr(cr i crpCr i 'CT^')Tr(cr :/ CT Al 'Cr : ,'Cr^). (A8) 



Finally, as a consequence of time-reversal symmetry, 
<r a G^(— k)<7 a = Goj(k) which implies X^ = Xc and 
Dr(w) = Gr(w).— Therefore, the cooperon and diffu- 
son contributions are equal and Eq. (|T5|) is obtained. 



Appendix B: Equal Rashba and Dresselhaus 
couplings 



The result for this case was discussed in Sec. IIVCI in 
terms of a simple unitary transformation. We show here 
that Eq. (J15I) is in agreement with this exact argument, 
even if Xo(q, u>) is approximated as in Eq. (fT8|) . By 
using 17 1 k) = ^X)±(l ^ cr y)l' c T ma + e x ) we can write 

GUk) = ^(l± ff „)Gi°'(kT)na + e s ), (Bl) 

z ± 

with GL (k) the Green's function in the absence of spin- 
orbit coupling. Inserting this formula in Eq. (|17[) it is 
easily obtained that, at a- = 0, 

+Tr[cr M (l ± <7„)<7„(1 T ^)]4 0) (qT2ma + e IjW )}. (B2) 

The value of XijXi'j' without spin-orbit coupling of 
Eq. (Jig was obtained by using A{, 0) (q,w) ~ l-(q£) 2 /2- 
|w|r. Substituting this expression in (|B2|) above, we find 
the same result of Eq. {JTSJ) at a_ = 0. Therefore, Eq. (|Tg|) 



satisfies the unitary transformation argument and direct 
numerical evaluation of Eq. (|15l) has to be in agreement 
with the discussion of Sec. IIVCI 



Appendix C: Details on the evaluation of XijXi'j' 
with only Rashba coupling 

To evaluate XijXi'j' it is convenient to rescale in 
Eq. (fl~5|) the frequency and the diffuson correlator as 



io = 2ljt(R/£) 2 , d R {w) = 2irmT£ 2 D R (Lj). (CI) 

such that the tensor amplitude appearing in Eq. (1251) is 
given by 

Vij.i'j' = 35 / wdw Yl d R v {w)d R " (w) 

fJb,V^' ,v' 

x r Tr(a i a l j i aiiaf J ,>)Tr(aja 1/ >aj>a v ). (C2) 

It is also convenient to use the dimensionlcss variable 
k = qi? in the Fourier transform (fTB]) 



ndndfj) 

tt£ 2 



-IK COs(<f> — (p) 



x #T t 1 ~ X d ( K /R> w{1/R) 2 /2t)] , (C3) 



dn(w) 



where it is found that the second line of Eq. (|C3[) only 
depends on the two ratios R/X±. Therefore, also f]ij,i'j' 
only depends on these two parameters, beside ip (the az- 
imuthal angle of R) . For example, in the case of Rashba 
spin-orbit coupling ((3 — 0), we can assume p = and 
perform the angular integration in terms of Bessel func- 
tions Ji(n) to obtain 

/d° K °(w) 

d»(w)- HndA ° dt{w) ° d " H 

[) ~J ° ° <CM 

V -d% z (w) d z K z (w), 



with 



(C4) 



<CH 



CM 



Jq{k) 

K 2 + W 



df{w) = 



4r k J\ (k) 
f(r, k, w) 



{k 2 + W + 4r 2 )Jg(K) 

f(r,K,w) 



d ±, w) = [f(r, k, w) + Sr 2 K 2 } Jq(k) T 8r 2 K 2 J 2 { K ) 
f(r,K,w)(n 2 +w + Ar 2 ) 

where r = R/ X and 

/(r, k, w) = (k, 2 + w) 2 - Ar 2 {K 2 - 3w) + 32r 4 . (C5) 

To obtain r\ij,i>y from Eq. (|C2[) . the remaining integra- 
tions (in dn and dw) are evaluated numerically We 
note that not all r\ij.i'y are independent. Some gen- 
eral symmetry transformations of rjij,i>j> are discussed 
in Appendix [D] but additional equalities exist in the 
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case of Rashba spin-orbit coupling, due to the simpli- 
fied structure of the diffuson matrix (|C4|) . For exam- 
ple, r] xx ,x x = Vzz, Z z and r) xx>xz = r) xz>xx , as discussed in 
Sec. IV CI Furthermore, several components are zero (in 
particular, r\i y ^y = if j' = x, z). The independent com- 
ponents of r/ijjiji are obtained in terms of the integrals 

Jo ^r^r." w dw, for example Vyy,yy = ^o is given by 



Vyy,: 



3 

8. /(i 



E[4M] 2 + 2ffiWf ^. 



(C6) 
where the sum runs over i = 0,x,y,z. 

Due to their special relevance, we give below the ex- 
pressions following from Eqs. (|38l) . (|44|) . and (J49J) for <fr, 
$' and $ ± : 



tan 2$ 



tan 4$' 



tan 2$+= 



tan 2$~ = 



2fZ°d%(d%-d!g)wdw 

±I™d%{d% + dU)wdw 
/o 00 W + ^) a -4(d^)>d«; , 

2f™d™d%wdw 
f °°d™(d% + d%)wdw' 

2 r o °° d^d^wdw 



(C7) 



(C8) 



(C9) 



(CIO) 



As seen from Eq. ()C4|) , the diffuson matrix can be diago- 
nalized at each value of w with a rotation around e y but 
the angle is determined by the ratio d^ /(dS" + d^) and 
in general depends on w. However, the results of Sec. [V] 
imply that such ratio is approximately independent of w. 
By writing 



2d n 



"R 



tan<f>. 



(Cll) 



where $ is a function of R/X (but not of w), Eqs. (IC7|) - 
(|C10|) can be simplified to give $ ~ $' ~ $ ± . 

As a final remark, we note that Eq. (|C11| also implies 
Aq ~ Ag and therefore Eq. (I5T1) for the variance. This 
is because the matrix dn(w) becomes approximately di- 
agonal with new rotated frames for Ij 2- By denoting as 
dn(w) the transformed (diagonal) diffuson, we have 

q r°° _ 

4) = | y o [«°) 2 + (d%) 2 + m 2 + (d%) 2 ]wdw, 

^ r°° _ 

^ = | y o [« ) 2 + to 2 - «) 2 - («) 2 m w , 

3 Z" 00 - - 
B * = l] dgdgwdw, 

3 f°° ~ - 

while the other coefficients of Eq. (f5Tj) are obtained from 
cyclic permutations of xyz. In particular, Aq = A' a 
within this approximation. 

Appendix D: Symmetries of »7y »/}•' 

One relation that the Tjijj'jt have to satisfy is simply 
due to XijXi'j' = Xi'j'Xij, which gives 



Vij,i'j' v / 'li , j , ,ij \ /* 



(Dl) 



A second relation can be obtained from relabeling of the 
magnetic moments 1 and 2 



7fa,i'j'(R) =Vji,j'i'(~'R)- 



(D2) 



Furthermore, the electron Hamiltonian is invariant upon 
a 7r rotation around z and we can transform the tensor 
r]jij'i'(— R) in Eq. (|D2j) to obtain 

?7y,i'i' (R) = Vjij'i' (TCjTZuTZjjTZi'i'Tlj'j' , (D3) 

where 72^ = fetf^ - ^x^jaj - 5 ly 5j V . 
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